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Abstract 

Background: During elongation, nnulti-subunit RNA polynnerases (RNAPs) cycle between phosphodiester bond 
fornnation and nucleic acid translocation. In the confornnation associated with catalysis, the mobile "trigger loop" of 
the catalytic subunit closes on the nucleoside triphosphate (NTP) substrate. Closing of the trigger loop is expected 
to exclude water from the active site, and dehydration may contribute to catalysis and fidelity. In the absence of a 
NTP substrate in the active site, the trigger loop opens, which may enable translocation. Another notable structural 
element of the RNAP catalytic center is the "bridge helix" that separates the active site from downstream DNA. The 
bridge helix may participate in translocation by bending against the RNA/DNA hybrid to induce RNAP forward 
movement and to vacate the active site for the next NTP loading. The transition between catalytic and 
translocation conformations of RNAP is not evident from static crystallographic snapshots in which macromolecular 
motions may be restrained by crystal packing. 

Results: All atom molecular dynamics simulations of Thermus thermophilus (Tt) RNAP reveal flexible hinges, located 
within the two helices at the base of the trigger loop, and two glycine hinges clustered near the N-terminal end of 
the bridge helix. As simulation progresses, these hinges adopt distinct conformations in the closed and open 
trigger loop structures. A number of residues (described as "switch" residues) trade atomic contacts (ion pairs or 
hydrogen bonds) in response to changes in hinge orientation. In vivo phenotypes and in vitro activities rendered 
by mutations in the hinge and switch residues in Socchoromyces cerevisioe (Sc) RNAP II support the importance of 
conformational changes predicted from simulations in catalysis and translocation. During simulation, the elongation 
complex with an open trigger loop spontaneously translocates forward relative to the elongation complex with a 
closed trigger loop. 

Conclusions: Switching between catalytic and translocating RNAP forms involves closing and opening of the 
trigger loop and long-range conformational changes in the atomic contacts of amino acid side chains, some 
located at a considerable distance from the trigger loop and active site. Trigger loop closing appears to support 
chemistry and the fidelity of RNA synthesis. Trigger loop opening and limited bridge helix bending appears to 
promote forward nucleic acid translocation. 
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Background 

In each phosphodiester bond addition cycle during RNA 
elongation, there are multiple intermediates, but two 
major events are catalysis (phosphodiester bond forma- 
tion) and translocation (movement of nucleic acids). A 
major question in transcription, therefore, is how RNA 
polymerase (RNAP) conformations and dynamics sup- 
port cycling between catalysis and translocation. From 
structures, it appears that a closed conformation of the 
"trigger loop" supports catalysis [1,2]. Two models have 
been put forward for translocation. In one model, in- 
coming NTP substrates act as allosteric effectors driving 
forward translocation [3-7]. In another model, thermally 
driven oscillations of RNAP elements stimulate nucleic 
acids to slide forward [8-11]. In this work, translocation 
is modeled and analyzed mostly in the presence of a 
NTP substrate, which may influence nucleic acid move- 
ments through RNAP. 

Two defining features of the multi-subunit RNAP tern- 
ary elongation complex (TEC) are the "bridge helix" and 
the "trigger loop", which are implicated in translocation, 
nucleotide sequestration and fidelity of NTP incorpor- 
ation [2,3,9,11]. The bridge helix, which abuts the active 
site and contacts the leading edge of the RNA/DNA 
hybrid, may function as a ratchet to drive forward trans- 
location, for instance, by bending against the RNA/DNA 
hybrid [12]. Thermus thermophilus (Tt) RNAP bridge 
helix (|3' 1067 to 1104) may have a glycine hinge (or two 
neighboring hinges) within the segment 1076-GARKGG- 
1081 [2,13-15]. Recent high throughput bridge helix 
mutagenesis of archaeal Methanocaldococcus jannashii 
(Mj) RNAP combined with molecular dynamics simula- 
tion supports the presence of two major bend points, cen- 
tered at Mj RNAP 809-GG-810 (corresponding to Tt 
RNAP 1076-GA-1077) and G825 (Tt RNAP G1092) 
[14,15]. In a Tt RNAP TEC with a closed trigger loop, mo- 
lecular dynamics simulations also support a primary bend 
near 1080-GG-1081 [13]. 

In the active site of the Tt RNAP TEC, the flexible trigger 
loop (|3' 1219 to 1265) adopts closed and open conforma- 
tions [1,16]. The RNAP TEC with a closed trigger loop has 
been observed in the presence of an accurately loaded NTP 
substrate and is thought to be the catalytic conformation. 
The open trigger loop TEC, thought to be the translocating 
form [17,18], was observed in the absence of a NTP 
substrate or in the presence of the antibiotic streptolydigin, 
which may inhibit transcription by allosterically blocking 
trigger loop closing, NTP substrate alignment or both. The 
RNAP trigger loop may have multiple roles in fidelity: 
maintaining the accuracy of RNA synthesis and maintain- 
ing the precision of the single base step length. Addition- 
ally, closing the trigger loop may dehydrate the RNAP 
catalytic center contributing to the fidelity of NTP recogni- 
tion [13]. Within a-helical segments surrounding the 



trigger loop, Tt RNAP has been suggested to have two 
GXP hinges [2,14,19], 1230-GEP-1232 and 1255-GLP-1257 
[20]. It is not clear whether these hinges are involved in 
translocation and/or phosphodiester bond formation or 
how motion of these hinges might be coupled to more dis- 
tant conformational changes of RNAP. 

We addressed the effects of trigger loop closing and 
opening on the RNAP TEC conformation using all atom 
molecular dynamics simulations, starting from known 
crystal structures (PDB files 2PPB, 205I and 205J) 
[16,21]. Molecular dynamics is a method to refine x-ray 
crystal structures in the presence of water and counter- 
ions and in the absence of crystal packing forces [22-24] . 
Additionally, molecular dynamics can explore otherwise 
unknown TEC conformations, and probe potential inter- 
actions between atoms. Recently, aU atom molecular dy- 
namics simulations of Saccharomyces cerevisiae (Sc) 
RNAP II were used to investigate trigger loop closing, 
NTP binding and pyrophosphate release [25,26]. For the 
current studies, Tt RNAP TECs were chosen over Sc 
RNAP II TECs for molecular dynamics analyses because 
of the higher resolution, better catalytic geometry and 
enhanced active site closing in the 205J Tt RNAP closed 
trigger loop TEC x-ray structure [16] compared to the 
2E2H Sc RNAP II structure [1]. Tighter active site clos- 
ing is expected to enhance the contrast between the 
closed and open Tt RNAP TECs in simulations. 

Modeling of Tt RNAP TECs with accurately loaded 
ATP substrates established that, as molecular dynamics 
simulations progress, closed and open TECs diverge into 
two distinct conformations that include bending of the 
bridge helix and trigger loop hinges and switching of 
multiple atomic contacts between amino acid residues 
that can be distal from the trigger loop and from the ac- 
tive site. The functional importance of hinge movements 
and the long-range conformational changes identified in 
the simulations were supported by site-directed muta- 
genesis and in vitro characterization of mutant RNAPs. 
Conformational changes during cycling between catalytic 
and translocating forms of RNAP indicate that closing 
and opening of the trigger loop and limited bridge helix 
bending may check and maintain the accuracy of trans- 
location steps. 

Results and discussion 

Molecular dynamics simulations of closed and open RNAP 
TECs 

We compare short duration molecular dynamics simula- 
tions of Tt RNAP with closed and open trigger loop 
conformations to gain potential insight into mechanism 
and to generate predictions for site-directed mutagen- 
esis. The RNAP TEC with a NTP bound and a closed 
trigger loop is thought to be the catalytic conformation, 
and the TEC with an open trigger loop is thought to 
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approach the translocating form. Simulations have been 
run several times with consistent results, so representa- 
tive trajectories are shown. A simulation of the closed 
catalytic TEC was done with restraints placed on the 3'- 
O of the RNA, the ATP substrate and Mg-I and Mg-II 
(referred to as "closed-Mg"). Restraining the active site 
to enforce catalytic geometry provides the greatest con- 
trast to the open TEC simulation and, therefore, helps to 
generate useful predictions. Simulations indicate signifi- 
cant conformational flexibility in RNAP without large 
structural distortions. In particular, modeling indicates 
distinct bending modes in the H1-H2 trigger loop and 
H3-H4 bridge helix hinges, supporting our description 
of these motifs (Figure 1). Closed and open TEC simula- 
tions diverge in conformation and atomic contacts, sug- 
gesting targets for RNAP mutagenesis of hinges and of 
residues nearby and at a distance that may support hinge 
movement. We find that the simulations are highly pre- 
dictive for the outcome of site-directed amino acid 
substitutions. 



Although strongly predictive for mutagenesis, current 
models have limitations. The RNAP mechanism is com- 
plex involving phosphodiester bond synthesis (closed 
trigger loop), pyrophosphate release (closed^open trig- 
ger loop), translocation (open trigger loop), NTP-Mg 
loading (open^closed trigger loop). Because the bridge 
helix borders the active site, helix bending and dynamics 
are likely to influence multiple steps, including catalysis 
and translocation, and also to influence off pathway 
steps, such as stalling, pausing, arrest and termination. 
The current simulations relate most closely to two 
mechanistic steps, phosphodiester bond synthesis and 
translocation. Current models are not adequate in their 
design or duration to address the energetics or kinetics 
of transitions between steps. The entire RNAP bond 
addition cycle occurs on a time scale of -1-50 millise- 
conds, so the current simulations (10.5 ns) provide lim- 
ited snapshots of stages of the mechanism. The models 
give insight into partial trigger loop closing and opening, 
active site hydration/dehydration during catalysis and 
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Figure 1 A model for trigger loop hinges, bridge helix hinges and bridge helix bending modes based on all atom molecular dynamics 
simulations. At the top of the figure, diagrams of the closed TEC, the closed product TEC (after chemistry) and the translocating TEC are shown. 
DNA is grey; RNA is red; the NTP substrate (or incorporated NMP and pyrophosphate) is blue; the trigger loop (TL) is purple; the bridge helix (BH) 
is yellow. Interpretations of simulations are shown schematically below. Simulations indicate trigger loop hinges HI and H2, bridge helix hinges 
H3 and H4 and bridge helix bend modes Bl (straighter) and B2 (more sharply bent). 
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thermal ratchet translocation, which appears to be 
strongly restrained by tight trigger loop closing during 
catalysis and enabled by trigger loop opening and bridge 
helix bending at an N- terminal hinge (s) (H3/H4). 

X-ray crystal structures of Tt RNAP TECs with a 
loaded ATP substrate and closed or open trigger loop 
conformations were analyzed by all atom molecular dy- 
namics simulations in explicit solvent including water 
and counterions. For simulations, AMPcPP (a-, p- 
methylene ATP) in the original crystal structures was 
substituted with ATP. The extent of convergence of 
simulations from RMSD (root mean square deviation 
compared to the original crystal structure) is indicated 
in Additional file 1: Figure SI. Additional file 1: Figure 
S2 shows that simulations maintain reasonable geom- 
etries of nucleic acids and stable DNA/DNA and RNA/ 
DNA base pairing. Most significantly, starting from very 
similar initial states, Tt RNAP TECs with closed or open 
trigger loops consistently diverged into distinct and 
recognizable conformations with characteristic altera- 
tions in functional atomic contacts. 

Similarly to DNA polymerases (DNAPs), RNAPs 
utilize a two Mg^"^ elongation mechanism [27,28]. In Tt 
RNAP, the metals at the active site are designated Mg-I 
and Mg-II, and in Sc RNAP II, the metals are designated 
Mg-A and Mg-B. It is likely that trigger loop closing, 
orientation of nucleotide substrates and associated 
metals, and dehydration of the active site are important 
aspects of DNAP and RNAP mechanism and fidelity 
[13,29]. A detailed model for the AMPcPP-Mg complex 
is presented in reference [16], which was used in our set 
up for a catalytic RNAP simulation with ATP substituted 
for AMPcPP. We find that binding Mg atoms to, 1) the 
RNA 3'-0, 2) ATP phosphate oxygens, and 3) key aspar- 
tic acid residues, results in a model that maintains cata- 
lytic geometry through simulation, as indicated by the 
RNA 3'-0 to NTP a-phosphate distance and the separ- 
ation of the active site Mg-I and Mg-II (Additional file 1: 
Figure S2B). 

In this work, three representative RNAP simulations 
are compared: 1) a closed trigger loop TEC with bonded 
Mg atoms (with restraints) and loaded ATP (referred to 
as "closed-Mg"); 2) a closed trigger loop TEC lacking 
bonded Mg atoms (without restraints) and with loaded 
ATP (referred to as "closed"); and 3) an open trigger 
loop TEC without bonded Mg atoms (without restraints) 
and with loaded ATP (referred to as "open"). The com- 
parisons, therefore, indicate the effects of restraining the 
catalytic center by bonding Mg atoms and also of trigger 
loop closing and opening. In the closed-Mg simulation, 
bonding of the Mg atoms was done to enforce a more 
precise catalytic geometry and also to maintain tighter 
trigger loop closing during simulation (Additional file 1: 
Figures SI and S2). For balanced comparison, a closed 



TEC simulation is shown without bonded Mg atoms. 
The closed TEC simulation maintains a closed trigger 
loop conformation but diverges from catalytic geometry 
(Additional file 1: Figure S2). In Additional file 1: Figure 
SI, RMSD analyses indicate that, as expected, the 
closed-Mg TEC simulation remains most similar to the 
original 205J closed trigger loop TEC crystal structure. 
Additionally, bonding the Mg atoms restrains the 
motions of nucleic acids and ATP substrate (Additional 
file 1: Figure SI). Motions of the nucleic acids are most 
similar for the closed-Mg and closed RNAP TECs (Add- 
itional file 1: Figures SI D-F), while overall motions of 
protein appear more similar for the closed and open 
trigger loop TECs (Additional file 1: Figures SID-F), 
demonstrating that both restraining the active site Mg 
atoms and closing of the trigger loop are consequential 
in simulation. 

A major conclusion of the molecular dynamics simula- 
tions is that the closed-Mg RNAP TEC shows significant 
conformational differences from the open TEC. Closed 
(205J) and open (2PPB) Tt RNAP TEC starting struc- 
tures differ somewhat in the position of the AMPcPP 
substrate analogue (AMPcPP is found in insertion 
(closed; 205J) and pre-insertion (open; 2PPB) sites), the 
downstream DNA clamp, the "lid" and "rudder". Overall, 
however, the original closed and open trigger loop TECs 
have a relatively straight conformation of the bridge 
helix and very similar conformations overall [16]. Not- 
ably, as simulations progress, much more dramatic dif- 
ferences develop than are present initially, comparing 
closed and open trigger loop RNAP TECs (Figure 1). 
Significant differences between these two conformations 
can be visualized readily in Movies Ml and M2. Sup- 
porting these distinct conformations, multiple ion pairs 
and hydrogen bonds in the two largest RNAP subunits 
switch atomic contacts in response to the conformation 
of the trigger loop. As expected, the closed TEC lacking 
bonded Mg atoms acquires an intermediate form that is 
generally most similar in conformation to the closed-Mg 
RNAP TEC. 

At the top of Figure 1, three TEC intermediates are 
shown schematically: 1) a catalytic TEC with NTP sub- 
strate loaded; 2) a product TEC after phosphodiester 
bond formation but before trigger loop opening; and 3) 
a translocating TEC with an open trigger loop and bent 
bridge helix. The closed-Mg simulation appears to be a 
reasonable representation of the catalytic TEC. The 
closed TEC deviates from catalytic geometry during 
simulation and may represent a marginally opened inter- 
mediate in bond addition (i.e. before phosphodiester 
bond formation). The open TEC simulation is close to 
the translocating TEC except that our model includes an 
accurately loaded ATP substrate, which may restrict for- 
ward RNA/DNA translocation. After simulation, the 
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closed-Mg TEC has a straight conformation of the 
bridge helix (Bl mode) with stretching between the H3 
and H4 hinges, a tightly closed trigger loop and reason- 
able catalytic geometry (Additional file 1: Figure S2). 
The open TEC has a more bent conformation of the 
bridge helix (B2 mode), a relaxed and open trigger loop 
and a tendency for forward translocation (compare 
Movies Ml and M2). 

RNAP translocation 

Previously, normal mode analysis and molecular dynamics 
simulations of Sc RNAP II TECs with closed and open trig- 
ger loop conformations were compared to obtain insight 
into RNAP translocation [17,18]. Because the structure of a 
pre-translocated Tt RNAP TEC with an open trigger loop 
has not yet been reported, a post-translocated Tt RNAP 
TEC was selected with an open trigger loop and NTP in the 
active site (2PPB). The open RNAP TEC was considered to 
be a reasonable comparison to the closed-Mg and closed 
TEC simulations (205J), which also have a NTP loaded In 
Figure 2, Additional file 1: S3 and S4, we compare the spon- 
taneous translocation of nucleic acids during simulation of 
the closed-Mg, closed and open trigger loop Tt RNAP TECs. 



To determine the extent of nucleic acid displacement in each 
TEC, translocation vectors were defined for downstream 
DNA/DNA and upstream RNA/DNA (Figure 2A) [18]. Sig- 
nificantly, the downstream DNA/DNA is translocated for- 
ward in the open TEC compared to the closed-Mg TEC. As 
a dramatic example, the i + 2 downstream base-pair appears 
to be a full step forward (>3 A) in the open TEC relative to 
the closed-Mg TEC. Notably, the upstream RNA/DNA re- 
gion is relatively compressed in the closed-Mg TEC and 
more extended in the open TEC. In Additional file 1: Figure 
S3, relevant overlay snapshot comparisons of closed-Mg, 
closed and open TEC simulations are shown to indicate 
translocation dynamics during simulation. 

In Figure 2B, a scatter plot is shown to compare trans- 
location positions of closed-Mg, closed and open RNAP 
TECs. Consistent with snapshots (Figure 2A and Add- 
itional file 1: Figure S3), the closed-Mg TEC is the least 
forward translocated for downstream DNA/DNA (~ -2 A 

on average) and upstream RNA/DNA ( 0.5 A) relative 

to initial crystal structures. Significantly, the open TEC is 
forward translocated for both downstream DNA/DNA 
(-- 3 A) and upstream RNA/DNA (-- 2 A). The closed TEC 
is intermediate between the closed-Mg TEC and the open 
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Figure 2 The open trigger loop RNAP TEC appears to be tiie translocating conformation. Vectors describe downstream DNA/DNA and 
upstream RNA/DNA translocation during simulation. The downstream translocation vector is drawn through DG5 (DNA template dGMP 5) and 
G37 (RNA GMP 37) 3' carbons. Translocation is measured as the displacement of an orthogonal projection of DA40 (non-template DNA dAMP 40) 
3' carbon onto the vector relative to the initial structure. The upstream translocation vector is drawn through DTI 4 (DNA template dTMP 14) and 
C31 (RNA CMP 31) 3' carbons. Translocation is measured as the displacement of an orthogonal projection of G33 (RNA GMP 33) 3' carbon onto 
the vector relative to the initial structure. Closed-Mg (green) and open (blue) TECs are shown. B) Scatter plot of upstream and downstream 
translocation (each point represents a 20 ps interval). "X" indicates 10.5 ns time points for the three simulations relative to a (0,0) starting position 
in the original crystal structures (205J for closed; 2PPB for open). 
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TEC in downstream DNA/DNA translocation (-OA) and 
similar to the open TEC in upstream RNA/DNA expan- 
sion (~ 2.5 A), emphasizing the role of the closed trigger 
loop in restricting translocation of the downstream DNA. 
Additional file 1: Figure S4 shows the trajectory of trans- 
location vectors as a function of simulation time. Closing 
of the trigger loop has the effect during simulation of 
restricting the forward translocation of the downstream 
DNA/DNA duplex. Restraining the active site in the 
closed-Mg simulation has the effect of restricting the 
stretching of the upstream RNA/DNA hybrid. 

The RMSD analysis shown in Additional file 1: Figure 
SI also indicates that nucleic acids are much more mo- 
bile in the open trigger loop RNAP TEC than in closed 
TECs, consistent with the open trigger loop enabling 
forward translocation. These results are in excellent 
agreement with those of Feig and Burton for Sc RNAP II 
[17,18], although a closed-Mg TEC with active site 
restraints, which shows the greatest contrast from the 
open TEC in this work, was not simulated in the previ- 
ous work. Taken together, our findings confirm the 
greater tendency for DNA/DNA and RNA/DNA forward 
translocation in the open TEC, which was predicted 
based on the previous molecular dynamics and crystal 
structure analyses [18]. In the open RNAP TEC, the 
RNA/DNA hybrid may be restrained from taking a full 
forward translocation step because of the presence of 
ATP, rather than incorporated 3 'AMP and byproduct 
pyrophosphate, in the active site. In the open TEC simu- 
lation, the bridge helix sharply bends in a B2 conform- 
ation (Figure 1; Movie M2) and holds a B2 bend without 
oscillation. Simulation, therefore, may indicate a bridge 
helix bend associated with at least partial forward trans- 
location. A similar conformational change has been sug- 
gested to initiate forward translocation of bacteriophage 
T7 RNAP after pyrophosphate release [30]. 

RNAP hinges 

Based on the flexibility of glycines and the helix breaking 
tendency of prolines, we propose that helices surround- 
ing the flexible trigger loop include hinges HI and H2 
(Figures 1 and 3). HI is predicted to lie within the se- 
quence |3' 1230-GEPGTQ-1235, and H2 is predicted 
within the sequence |3' 1255-GLP-1257 [31]. The simula- 
tion reveals bending and hinge dynamics in the vicinity 
of G1230 (HI) and G1255 (H2). Movies Ml and M2 
show trigger loop, bridge helix and nucleic acid dynam- 
ics for the closed-Mg and open RNAP TEC simulations. 
The movies provide a means to visualize motions of the 
HI and H2 hinges. In Supplementary Materials, we show 
trajectory data (simulation time versus a dynamics 
metric (either Apsi angle (a measure of amino acid resi- 
due conformation) or secondary structure)) indicating 
dynamic actions of HI (Additional file 1: Figure S5) and 



H2 (Additional file 1: Figure S6). These graphs and the 
movies indicate that HI and H2 are complex hinges with 
multiple modes of bending and dynamics depending on 
the state of the RNAP TEC and its position in the bond 
addition cycle. 

The trigger loop is in close association with the bridge 
helix [9]. Molecular dynamics simulations indicate that 
Tt RNAP p' 1076-GARKGG-1081 on the bridge helix 
contains two nearby and semi-independent hinges 
described here as H3 (1076-GA-1077) and H4 (1080- 
GG- 1081) (Figure 1). Flexibility in and around H3 and 
H4 glycines results in distinct modes of bridge helix 
bending referred to as Bl (closed-Mg TEC; a relatively 
straight bridge helix conformation) and B2 (open TEC; a 
strongly bent conformation) (Figures 1 and 4). The 
strongly bent B2 mode in the open TEC in particular 
might be expected to promote forward translocation of 
RNAP, consistent with an open TEC being a translocat- 
ing conformation [17,18] (Figure 2, Additional file 1: S3 
and S4). Distinct bridge helix bending modes Bl and B2 
develop early in simulation and are stably maintained 
(Additional file 2: Movies Ml and Additional file 3: M2). 
In Additional file 1: Figure S7, trajectories of Apsi versus 
simulation time and of secondary structure versus simu- 
lation time are shown for bridge helix hinges H3/H4. 
Additional metrics that indicate distinct bending and dy- 
namics of trigger loop hinge HI and bridge helix hinges 
H3/H4 comparing closed and open RNAP TECs are 
shown in Additional file 1: Figure S8. Additional file 1: 
Figure S8A indicates different bending modes of the HI 
hinge comparing the closed-Mg to the closed and open 
TECs. Additional file 1: Figure S8B indicates different 
bending modes of the H3/H4 bridge helix hinges in the 
three simulations. We propose that motions and bend- 
ing of hinges HI to H4 initiates a larger conformational 
switch between catalytic (closed-Mg) and translocating 
(open) forms during each RNAP phosphodiester bond 
addition and nucleic acid translocation cycle. Slight 
bending of HI and H2 coincides with acquisition of the 
catalytic conformation by the closed-Mg TEC. Mobility 
of the H2 hinge is pronounced both in the catalytic and 
translocating conformations, but may be specifically im- 
portant for trigger loop opening (note the loss of the 
helical structure starting from H2 in the open versus 
closed conformations in Figure 3) [16,21]. 

RNAP hinges in transcription 

Because molecular dynamics simulation predicted that HI 
to H4 are dynamic hinges, conserved hinge residues in 
Saccharomyces cerevisiae (Sc) RNAP II were substituted 
by in vitro mutagenesis and assayed for transcription rate 
and fidelity (Figure 5A). Mutations were made in Sc 
RNAP II rather than Tt RNAP for both pragmatic and 
tactical reasons. From a practical standpoint, we were 
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Closed-Mg 



1219-EAVGIVAAQSIGEPGTQLTMRTFHTGGVAGAADITQGLPRVIELFE-1264 Tt b' 
913-EAIGVIAAQSIGEPGTQLTMRTFHIGGAASRA-ITGGLPRVADLFE-1146 Ec b' 
5 8-EAVGIVAAQSIGEPGTQMTMRTFHYAGVAEIN-VTLGLPRMIEIVD-5 73 Mj A" 
10 62-EMVGVLAAQSIGEPATQMTLNTFHFAGVASKK-VTSGVPRLKEILN-110 6 Sc Rpbl 

Trigger loop 



H1 



Figure 3 The trigger loop hinges HI (1230-GEPGTQ-1235) and H2 (1255-GLP-1257). Snapshots of the trigger loop region in open, closed 
and closed-Mg TECs are compared at about 10 ns simulation. Coloring is in "secondary structure" representation: helix is magenta, 3io helix is 
blue, coil is white, turn is cyan, (3-sheet is yellow. G1230, PI 232, HI 242, G1255 and PI 257 are in stick representation to locate the HI and H2 
hinges and orient the view. ATP and Mg are in space filling representation to locate the active site. Blue arrows indicate differences in secondary 
structure that develop during simulation. 



,D429 




H3/H4 



10 67-VLEYFISSHGARKGGADTALRTADSGYLTRKLVDVTHE-1104 Tt b' 
769-VLQYFISTHGARKGLADTALKTANSGYLTRRLVDVACP-806 Ec b' 
800-PTEFFFHA]yiGGREGLVDQAVRTAQSGYMQRRLINALg)-837 Mj A' 
810-PQEFFFHA]yiGGREGLIDTAVKTAETGYIQRRLVKALED-84 7 Sc Rpbl 
Bridge helix 



42 0-RERAGFDVRDVHRTHYGR-43 7 Tt b 
54 0-RERAGFEVRDVHPTHYGR-55 7 Ec b 

43 6-RSQPHFEARELHGTHWGK-453 Mj B" 

5 04-RDGKLAKPRQLHNTHWGL-521 Sc Rpb2 
Fork 



Figure 4 The bridge a-helix hinges and bending. To the right, 10 ns snapshots of closed, closed-Mg and open bridge helices are compared, 
colored in secondary structure representation, as in Figure 3. G1076, G1080 and G1081 are in space filling representation so that H3 and H4 
hinges can be identified. To the left, proposed switch residues (3' HI 075, K1079, D1083, R1087, D1090 and (3 D429 are shown in stick 
representation for closed-Mg and open TEC bridge helix structures. Blue arrows indicate differences in secondary structure that develop during 
simulations. 
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Figure 5 Transcriptional activities of hinge HI to H4 amino acid substitutions. A) Elongation and fidelity assays. Blue bars represent 
transcriptional elongation rate relative to wild type RNAP II at 200 NTP concentration. Error bars indicate standard error from exponential 
curve fitting. Red bars represent transcriptional fidelity, determined by a competition assay (Additional file 1: Figures S9 and SIO). Error bars are 
standard deviation of at least three independent determinations. B) Conserved ionic interactions surrounding the Sc RNAP II H2 hinge. 



initially better equipped to make the substitutions in Sc 
RNAP IL Furthermore, testing residues that are conserved 
across species provides more general insight into RNAP 
catalytic and translocation functions. Transcription rates 
are compared to wild type Sc RNAP II based on quantifi- 
cation of bulk elongation assays. The transcriptional fidel- 
ity analysis is based on a competition assay, in which 
AMP from ATP (the incorrect substrate at a high concen- 
tration) is misincorporated for GMP from GTP (the ac- 
curately templated substrate at a limiting concentration) 
(Additional file 1: Figures S9 and SIO)) [32]. 

Mutations in the trigger loop hinge HI (Sc RNAP II 
Rpbl 1073-GEPGTQ-1078) were expected to mainly 
affect formation of the catalytic conformation of the 
TEC because of its close proximity to the substrate NTP 
binding site and distinct conformational changes of the 
closed-Mg TEC observed in molecular dynamics simula- 
tions. Indeed, the G 1073 A substitution is very slow in 
elongation and demonstrates a significant fidelity defect 
(Figure 5A). Remarkably, P1075A, T1077A and Q1078A 
[2,33,34] are lethal substitutions in yeast, indicating es- 
sential functions. In Tt RNAP, Q1235A (corresponding 
to Sc RNAP II Q1078) is very slow in elongation and 
shows a strong fidelity defect [29]. G1076 substitutions 
were not made in this study because all possible hom- 
ologous substitutions were previously reported for Mj 
RNAP and were shown to have a moderate effect on 
transcription [15]. Interestingly, E1074A appears faster 
than wild type Sc RNAP II in elongation, suggesting that 
E1074 might play a non-essential regulatory role that 
appears dispensable for maintenance of the catalytic 
conformation. The suggested regulatory role of El 074 
does not appear to be related to mechanisms of cognate 
NTP selection, because this mutation does not affect 



transcriptional fidelity. Taken together, the results of HI 
mutagenesis confirm the important role of the HI hinge 
region in phosphodiester bond formation. 

As mentioned above, hinge H2 appears to be important 
for trigger loop opening. Notably, delayed trigger loop 
opening has been proposed to impair NTP selection, de- 
creasing transcriptional fidelity [2,10,20,34]. In agreement 
with the previously reported structural and functional evi- 
dence, substitutions in and around trigger loop hinge H2 
of Sc RNAP II tend to be rapid in elongation and some 
show significant transcriptional fidelity defects. Sc RNAP 
II has the sequence 1097-GVPRLKE-1103, with major H2 
bending predicted near to G1097. Significantly, substitut- 
ing 1097-GVP-1099 to 1097-AVA-1099 generates a Sc 
RNAP II that is rapid in elongation with a moderate fidel- 
ity defect (Figure 5A). As previously reported and as con- 
firmed here, E1103G, located near H2 is strongly defective 
in transcriptional fidelity [2,10,20]. El 103 forms an ionic 
contact to RllOO, which forms an ionic contact to E1351 
(Figure 5B). Similarly to the 1097-AVA-1099 H2 double 
substitution, RllOOA, E1103G and E1351A appear to be 
rapid in elongation. RllOOA and E1351A substitutions, 
however, do not render RNAP II error-prone, suggesting 
that transcription fidelity regulation may be restricted 
within the H2 hinge region, which includes G1097, P1099, 
RllOO and El 103. The gain of function (fast) phenotype 
of all the substitutions in and around H2 tested here sug- 
gests that transcription elongation rate may be regulated 
not only by H2 but also perhaps by downstream DNA 
contacts, which may be affected by H2 mobility. Because 
both G1097 and P1099 can be replaced with alanine, the 
mobility and bending of the H2 hinge observed in mo- 
lecular dynamics simulations does not appear essential for 
formation of a catalytic RNAP conformation but rather 
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may regulate the rate of phosphodiester bond formation 
and appears to play an important role in transcriptional 
fidelity. 

All atom molecular dynamics simulations of Tt RNAP 
suggest that bridge helix hinges H3 and H4 (Sc RNAP II 
Rpbl 819-GGREGL-824) undergo extensive conformational 
changes in transition from catalysis to translocation. Sc 
RNAP II Rpbl G819A and G820A are slow in elongation 
with moderate fidelity defects (Figure 5A), indicating that 
mobility of H3 contributes to an appropriate catalytic con- 
formation. R821 is mounted between the clustered H3 and 
H4 hinges and connects the bridge helix with the Rpb2 fork 
region. R821A is a lethal substitution. R821Q is viable but 
slow in elongation. The phenotypes of R821 alleles suggest 
the functional importance of bridge helix interactions with 
the RNAP "fork" region (Tt RNAP p 420 to 450), which 
approaches the active site. G823A appears faster than wild 
type Sc RNAP II in elongation. The 819-GGREGL-824 ^ 
819-AAREAL-824 triple alanine for glycine substitution in 
H3/H4 is a lethal substitution, suggesting the importance of 
flexibility and of the ability to deviate from precise a-helical 
structure in the H3/H4 hinge region. Extensive mutagenesis 
of the bridge helix H3/H4 region has been reported for Mj 
RNAP [14] and alanine scanning has been done for Escheri- 
chia coli (Ec) RNAP [35]. These previous studies, which also 
reported both fast and slow elongation phenotypes ren- 
dered by mutations in bridge helix residues are very consist- 
ent with the results reported here for Sc RNAP II H3/H4. 

To test for possible synergy between RNAP trigger loop 
and bridge helix hinges, two additional substitutions were 
made. The Rpbl G819A, G1073A double substitution 
combines a mutation in H3 (bridge) with one in HI 
(trigger loop). In transcriptional activities, the effects of 
combining these two substitutions appear to be additive 
rather than synergistic. Similarly, the effects of combining 
G819A, G1097A, and P1099A substitutions in H3 (bridge) 
and H2 (trigger loop) also appear to be approximately 
additive in transcriptional effects. These results are most 
consistent with hinges on the trigger loop and the bridge 
helix being largely independent of one another in the 
RNAP bond addition cycle, because tight coupling of trig- 
ger loop and bridge helix hinge motions might be 
expected to produce synergistic effects of combining hinge 
mutations. In the work of others, multiple RNAP muta- 
tions have been combined around the bridge helix C- 
terminal glycine hinge of Mj RNAP (G825) [15] and 
within the trigger loop region of Sc RNAP II [2]. 

Exonuclease III mapping 

To further investigate possible interactions between 
RNAP hinges, exonuclease III mapping was done from 
the upstream RNAP border for a 3' chain-terminated 
All TEC in the absence and presence of the incoming 
NTP (Figure 6). A schematic representation of the 3 '^5' 



exonuclease III back border mapping experiment is 
shown at the bottom of Figure 6A. Mapping provides 
insight into the resting translocation states of the Sc 
RNAP II TEC and the effects of the incoming cognate 
NTP on translocation [20,36]. There is a strong correl- 
ation between the stable sequestration of the cognate 
NTP in the RNAP II active site and the magnitude of the 
incoming NTP effect on the pre- to post-translocated 
border transition in the exonuclease III assay [20,36], 
showing that the post-translocated state is stabilized by 
binding of the cognate NTP in the active site [8] . Alterna- 
tively, the apparent stimulation of translocation might re- 
flect an allosteric effect of the incoming NTP previously 
proposed for mammalian RNAP II [4,37] and Ec RNAP 
[6,7,38]. 

In the experiment shown in Figure 6, the cognate CTP 
(encoded at C12) promotes the pre-to-post-translocated 
border transition in 3'dAll chain terminated Sc RNAP II 
TECs (Figure 6A, compare lanes 1-3 and 4-6 and 
Figures 6B-C; note that in Figure 6C the approximately 2- 
fold increase of the post-translocated state for wild type 
RNAP II corresponds to the increase of the post- 
translocated fraction from 25% to 51% as shown in 
Figure 6B). G819A substitution targeting hinge H3 does 
not affect RNAP II interaction with the incoming CTP 
(Figures 6B-C). The G1073A substitution in the trigger 
loop hinge HI leads to an apparent decrease of the incom- 
ing CTP effect on translocation. This effect may be due to 
the post- translocated shift of 3'dAll TEC equilibrium as 
compared to the wild type TEC (Figure 6B, compare the 
wild type and G1073A RNAP II without CTP), but not to 
a defective interaction of G1073A RNAP II with the in- 
coming CTP (Figure 6B, compare the wild type and 
G1073A in the presence of CTP). Notably, the combin- 
ation of mutations in the trigger loop hinge HI and bridge 
helix hinge H3 practically abolishes the ability of the in- 
coming CTP to stabilize the post-translocated state 
(Figure 6A, compare lanes 4-6 and 10-12; also Figure 6B- 
C). Because these two mutations have apparently additive 
effects on the elongation rate and fidelity of RNAP II 
(Figure 5A), their combined effects on NTP stabilization 
of the post-translocated TEC (Figure 6) may reflect the 
complexity of the RNAP elongation cycle or allosteric 
effects of NTPs. The translocation assay in the presence of 
the cognate NTP indicates a synergistic action of HI and 
H3 hinges in NTP binding. Unlike HI, the distal trigger 
loop hinge H2 does not appear to have synergy with 
bridge helix hinge H3 in the exonuclease Ill-based trans- 
location assay. In agreement with the results of the elong- 
ation and fidelity assays, which demonstrate a gain of 
function phenotype in elongation and decreased transcrip- 
tional fidelity, mutations in trigger loop hinge H2 
(G1097A, P1099A) increase the effects of the incoming 
CTP on the translocation equilibrium, further supporting 



Kireeva et a I. BMC Biophysics 201 2, 5:1 1 
http://www.bionnedcentral.conn/2046-1 682/5/1 1 



Page 10 of 18 



Time with 
exo III, s 

DNA 

Post^ 

RNA 3'dAll 



WT 

CTP 

20 40 90 20 40 90 



Rpb1 -G81 9A Rpb1 -G1 097A 
G1073A P1099A 

CTP CTP 

20 40 90 20 40 90 20 40 90 20 40 90 



1 2 3 4 5 



7 8 9 10 11 12 13 14 15 16 17 18 



dC49 
dC48 



Upstream 



SAUCGAGAGGGA^ 

3 ' GGATATCCTATGAACC|^CGTAGCTCTCCC^TGCCGCTTATGGGTA( 



TDS 

JTTAGCCGGACGACCA5'* 



1^ 




Figure 6 Exonuclease III mapping of Sc RNAP II TECs. A) Back border (upstream border) 3'-dAl 1 Sc RNAP II TEC mapping. 3'-dAl 1 (3'-dATP 
chain terminated Al 1 TECs) RNA controls, run in the same reaction, are shown to demonstrate effective chain termination at the Al 1 position. 
The strategy for mapping Sc RNAP II TECs in the absence or presence of incoming CTP is indicated at the bottom of the figure. B) 
Phosphorimager quantification of triplicate determinations done at 30 s exonuclease III incubation. C) Percent increase in the post-translocated 
state induced with CTP compared to wild type Sc RNAP II, from the data shown in B. Error bars represent standard deviation. 



the importance of H2 in cognate NTP selection by RNAP 
[20]. The significant effect of the bridge helix hinge H3 
G819A substitution on translocation in the presence of ei- 
ther a trigger loop hinge HI G1073A single or a H2 1097- 
GVP-1099 ^ 1097-AVA-1099 double substitution 
(Figure 6B-C) indicates that the bridge helix H3 hinge re- 
gion may have a role in translocation, as indicated in 
simulations (Movies Ml and M2). 

Importantly, the effects of mutations in the hinges on 
translocation have been pronounced only in the transloca- 
tion assays performed in the presence of the incoming 
NTP. This result may indicate that mutations in the 
hinges largely affect trigger loop and bridge helix interac- 
tions with the NTP positioned for catalysis in the active 
site of the closed TEC. Alternatively, it cannot be excluded 



that translocation of the stalled TEC deprived of NTPs 
and translocation during ongoing transcription, which 
takes place in the presence of the incoming NTP, occur by 
distinct pathways. Notably, very recent optical trapping 
single molecule experiments suggest that the NTP can 
bind both the pre- and post-translocated TECs [10]. 

Hinge conformations appear to be supported by 
alternate switch residue contacts 

Different modes of hinge HI to H4 action and Bl and B2 
bridge helix bending result in long-range changes in the 
atomic contacts of amino acid residues of Tt RNAP p and 
|3' subunits. We refer to the RNAP amino acids making 
altered contacts as "switch" residues. Note that these resi- 
dues should not be confused with disordered segments of 
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RNAP chains that assume increased structural order upon 
nucleic acid binding (those segments of chain previously 
designated "switches 1-5") [12,39]. In this paper, switch 
residues are amino acid side chains predicted from mo- 
lecular dynamics simulation to alter their molecular con- 
tacts (ionic contacts or hydrogen bonds) in response to 
closing and opening of the trigger loop. 

Molecular dynamics simulation reveals a charge relay 
chain across the bridge helix involving Tt RNAP |3' 
H1075, K1079, D1083, R1087, D1090 and fork residue p 
D429 (Figure 4, left images). Additional file 1: Figures 
S11-S13 show atomic contacts and simulation trajectories 
that relate to the proposed switching mechanism. |3' 
K1079, which is mounted between the H3 and H4 hinges, 
and |3 D429, with which K1079 can interact (open TEC), 
appear to be central residues in supporting different con- 
formations of the bridge helix and interacting fork. Com- 
paring late simulation snapshots of closed-Mg and open 
TECs, K1079 switches its contacts between D1083 
(closed-Mg) and D429 (open TEC). Tt RNAP K1079 and 
D429 are not conserved in Sc RNAP II, indicating that the 
core of this switching mechanism is specific to bacterial 
RNAPs. Tt RNAP p' D1083, R1087 and D1090 in the 
charge relay chain, however, are conserved in evolution 
and are expected to fulfill a similar role in both Tt RNAP 
and Sc RNAP II. The homologous substitution to Tt 
RNAP p' K1079A has been made in Ec RNAP and is 



strongly defective in elongation [35], as we would predict 
based on the central role of this switch residue in simula- 
tions. R1078, adjacent to K1079 and mounted between the 
bridge helix H3/H4 hinges, is also important in regulating 
the structure of the fork (Additional file 1: Figures Sll- 
S13). Interestingly, Tt p' D1090 is substituted with aspara- 
gine in Ec RNAP and glutamine in Mj RNAP, changes that 
would diminish interactions with R1087. Notably, Mj 
RNAP Q^E/D substitutions at the position correspond- 
ing to Tt RNAP p' D1090 are strong gain of function 
mutants (rapid in elongation) [14], consistent with the 
predicted importance of Tt p' D1090 in supporting bridge 
helix bending and dynamics. 

Based on all atom simulations, predictions for the Tt 
RNAP fork region and its contacts can potentially be 
extended. Specifically, Tt RNAP p R142, R331, R420, 
D426 and D429 are predicted to make potentially signifi- 
cant switching contacts (Figure 7). R142, R331 and po- 
tentially R420 can interact with D426, which is close to 
D429, which in the open TEC appears to make a key 
interaction to p' K1079 on the bridge helix (Figures 7 
and Additional file 1: Figures S11-S12). We posit that 
maintaining different ionic contacts to D426 by R142, 
R331 and R420 reorients D429 to approach bridge helix 
residue H1075 (closed-Mg) or K1079 (open). If H1075 
can take on a +1 charge by protonation, H1075 could 
form an ion pair with D429 in the closed-Mg TEC. 



Closed-Mg 




136-IVSQIHRSP GVYF-148 Tt p 

145-IVSQLHRSP GVFF-157 Ec (3 

136-ELIDLGEDPEDPFGYF-151 Mj B" 
188-DLYKLKECPFDMGGYF-2 03 Sc Rpb2 
Protrusion 



323-DDIDHLGNRRIRTVGELMT-341 Tt P 

443-DDIDHLGNRRIRSVGEMAE-4 61 Ec p 

334-DDKDHYAYKRAKLAGDLME-3 52 Mj B" 

3 96-DDRDHFGKKRLDLAGPLLA-414 Sc Rpb2 



42 0-RERAGFDVRDVHRTHYGR-43 7 Tt p 

54 0-RERAGFEVRDVHRTHYGR-557 Ec P 

436-RSQPHFEARELHGTHWGK-453 Mj B" 

5 04-RDGKLAKPRQLHNTHWGL-521 Sc Rpb2 
Fork 



Figure 7 Proposed interactions between bridge helix switch residues P' HI 075 and K1079 and the foric region. Fork residues (3 R420, 
D426 and D429 are proposed from simulations to be switcli residues. Interacting residues (3 R142 and R331 also appear to be involved in the 
switching mechanism. 
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A conserved switch near the RNAP active site 

From simulations of open RNAP TECs, a conserved 
switch is predicted beneath the Tt RNAP active site near 
Mg-II centered about 20 A from the HI trigger loop hinge 
(Figure 8). Additional file 1: Figures S14-S17 show atomic 
contacts that may be altered comparing the closed-Mg 
and open TEC RNAP snapshots. Tt RNAP |3 E685 and 
D686 appear invariant in evolution and closely approach 
Mg-II in x-ray crystal structures [40-42]. The central resi- 
due in the predicted switch is |3' D784. We posit that a po- 
tential role of this switch may be in release of the RNAP 
reaction byproduct pyrophosphate (diphosphate) and Mg- 
II after phosphodiester bond formation. Most significantly, 
in all closed RNAP TEC simulations, D784 maintained a 
distinctly different set of molecular contacts (more similar 
to the initial x-ray crystal structure) than those that develop 
during simulation of the open TEC. The images in Figure 8 
are representative of multiple simulations. Remarkably, the 
carboxyl head group of f> D784 points up (in this view) to- 
ward the main chain NH of p E685 to form a hydrogen 
bond in the closed-Mg RNAP TEC. In the open TEC, by 
contrast, the carboxyl head group of |3' D784 is consistently 
switched downward toward the |3' S942 hydroxyl head 
group to form a hydrogen bond. These interactions involve 
contacts of |3' D784 with p N683 and p N872. Only in the 
open TEC, an ion pair develops between p D686 and p 



R557, after D686 and E685 dissociate from Mg-II (Figures 8 
and Additional file 1: S14-S17). These simulations, there- 
fore, predict that p' D784 and S942 and p R557, N683, 
E685, D686 and N872 are switch residues participating in 
interactions with Mg-II that alter atomic contacts between 
the catalytic and translocating conformations of the Tt 
RNAP TEC; p' D784 appears to be of central importance 
in operating the switch. 

Site-directed mutagenesis of homologous switch resi- 
dues identified in simulations and a number of nearby 
active site residues was done in Sc RNAP II (Figure 9). 
Consistent with the proposed central role for D784 in 
this switch, Sc RNAP II Rpbl D526A, corresponding to 
Tt RNAP p' D784, is a lethal substitution. Sc RNAP II 
Rpb2 N834A, corresponding to Tt RNAP p N683, is 
very slow in elongation, consistent with its proposed role 
supporting the switch. Sc RNAP II Rpb2 N1013A, corre- 
sponding to Tt RNAP switch residue p N872, is slow in 
elongation and appears to have enhanced transcriptional 
fidelity relative to wild type Sc RNAP II. 

Consistent with the results of others, Sc RNAP II Mg-B 
interacting residues Rpb2 E836 and D837 (corresponding 
to Tt RNAP p E685 and D686) substitutions have signifi- 
cant effects on transcription [40-42]. Rpb2 D837A is a le- 
thal substitution. Rpb2 E836A is slow in elongation and 
defective in transcriptional fidelity. The defect in E836A 




Switch 



770-LSPASGEPLAKPSRDIILGLY-790 Tt P' 
491-LSPANGEPIIVPSQDWLGLY-511 Ec p' 
949-LSPRFGGPIIGAIHDFISGAY-969 Mj A' 
512-VSPQSNKPCMGIVQDTLCGIR-532 Sc Rpbl 
Switch/ Funnel /Pore 



R557 

678-PFDGYNFEDAIVI-690 Tt P 
807-PWNGYNFEDSILV-819 Ec P 
217-SYEGYNMEDAIVF-229 Mj B' 
82 9-CYSGYNQEDSMIM 841 Sc Rpb2 
Switch/Mg-II binding 

545-NLIPFLEHEDANRALMGSNM-564 Tt P 
666-SLIPFLEHEDANRALMGANM-685 Ec p 
142-GVAPYPEHNSAPRITMAAAM-161 Mj B' 
754-SIIPFPDHNQSPRNTYQSAM-773 Sc 
Funnel/Active site Rpb2 



864 -GTPVDVILNPLGVPSRMNLGQILET-888 Tt p 
109 0-GTPVDIVLNPLGVPSRMNIGQI LET- 1115 Ec P 

393-GIVPDIIINPHAIPSRMTVGQLLEM-417 Mj B' 
1005-GIVPDLIINPHAIPSRMTVAHLIEC-1029 Sc Rpb2 
Switch/ Funnel /Pore 



Figure 8 A conserved conformational switching mechanism centered on (3' D784 about 20 A from HI. The substrate ATP, Mg-I and Mg-II 
(green) are in space-filling representation. Proposed switch residues are in stick representation. The two orientations of D784 observed in 
simulations of closed and open trigger loop RNAP TECs are indicated with blue arrows. Some atomic interactions are indicated with blue lines or 
emphasized with blue circles. 
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Figure 9 Transcriptional activities and fidelity of Sc RNAP II substitutions in the active site and proposed switch residues (as in Figure 5A). 



fidelity can be attributed to misalignment of Mg-B and, 
thus, the substrate NTP during phosphodiester bond for- 
mation. Other substitutions that are expected to misalign 
the NTP substrate are Rpb2 S1019A and S1019F, which 
are adjacent to R1020, which corresponds to Tt RNAP p 
R879 (Figure 8). Sc RNAP II Rpb2 R1020K/Q substitutions 
are lethal (Figure 9) [43], indicating an essential function. 
The bulkier S1019F substitution is more strongly defective 
in transcription and fidelity than the S 1019 A substitution, 
as expected, if adjacent R1020 orients the NTP substrate 
for catalysis through its contact to the NTP y-phosphate. 
Sc RNAP II Rpb2 R766A/Q substitutions, corresponding 
to Tt RNAP p R557, are lethal (Figure 9) consistent with 
simulations (Figures 8 and Additional file 1: S17) [25]. 

Other substitutions, addressing the functional roles of 
the Rpb2 subunit switch residues located in the funnel 
and pore regions of the Sc RNAP II TEC are described 
in more detail in reference [43]. R983 is expected to 
form an ionic contact with E1028, and G985 and G988 
ensure the flexibility of R983. Note that enhanced pro- 
pensity to undergo transcriptional stalling or arrest (i.e. 
by backtracking) observed in Rpb2 mutants R983A, 
G985A, G988A and E1028Q may indicate a defect in 
forward translocation of the stalled TEC. These muta- 
tions might also affect RNAP interactions with the 3' 
RNA extruded from the active site during backtracking, 
stabilizing the arrested conformation of the TEC. 

Active site dehydration in catalysis 

In x-ray crystal structures, often water that is present 
within a protein is not resolved. Molecular dynamics simu- 
lation in explicit solvent, however, generates a model for in- 
ternal hydration of a protein. Previously, trigger loop 
closing was posited to drive water from the RNAP active 
site and to order water to optimize conditions for catalysis 
[13]. This model for active site closing and dehydration to 



order water is potentially a powerful and general concept 
to describe catalysis within enclosed active sites of enzymes 
[44-47]. In support of such a model for RNAP active site 
dehydration during trigger loop closing, the closed-Mg, 
closed and open RNAP TECs are compared for active site 
hydration (Additional file 1: Figures S18-S20). At the ATP 
N3 position, the closed-Mg TEC remains strongly dehy- 
drated throughout the simulation (Additional file 1: Figure 
SI 8 A). The closed TEC, which deviates from catalytic 
geometry (Additional file 1: Figure S2) and may not be as 
tightly closed as the closed-Mg TEC, is slightly more 
hydrated. As expected, the open trigger loop, which admits 
more access to water, is much more strongly hydrated 
around ATP N3. Measuring water surrounding the RNA 
3'-0, similar results are obtained (Additional file 1: Figure 
S18B). We conclude that at least one function of trigger 
loop closing is to reduce the hydration of the NTP sub- 
strate and to enhance the ordering of waters in the active 
site. Tightly bound waters are detected at Mg-I in both 
closed-Mg (one bound water) and open (three bound 
waters) RNAP TEC simulations (Additional file 1: Figures 
S19-S20). 

Model for RNAP cycling between catalysis and 
translocation 

From all atom molecular dynamics simulations of Tt 
RNAP TECs, we posit that RNAP may cycle during 
elongation between at least two distinct conformations 
that support either catalysis or translocation. The two 
conformations are characterized by bending at hinges 
within helical segments surrounding the trigger loop 
(HI and H2) and toward the N- terminal end of the 
bridge helix (H3/H4) (Figures 1, 3 and 4). Hinge motions 
are supported by amino acid "switch" residues (distinct 
from RNAP "switch regions 1-5") that change atomic 
contacts in response to hinge action. Remarkably, these 
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switch residues can be located at significant distances 
from the hinges and from the RNAP active site. The 
catalytic TEC appears to resist forward translocation be- 
cause of the tight closing of the trigger loop over the 
NTP substrate. By contrast, the RNAP TEC with an open 
trigger loop appears to support increased forward nucleic 
acid translocation and enhanced nucleic acid motion. Fur- 
thermore, in the RNAP TEC with an open trigger loop 
conformation, the bridge helix becomes more strongly 
bent against the RNA/DNA hybrid in simulation, indicat- 
ing that mechanical pushing by the bridge helix may con- 
tribute to a fractional forward translocation step. In 
current simulations performed on a nanosecond time- 
scale, there is little indication of an oscillating RNAP 
ratchet that might drive both forward and reverse trans- 
location and unrestrained nucleic acid sliding. Rather, 
bridge helix bending against the RNA/DNA hybrid 
appears to generate a more deliberate but partial forward 
step. Very similarly, during movement of bacteriophage 
T7 RNAP, a steric clash develops between the pre- 
translocated RNA/DNA hybrid and the open product 
RNAP TEC stimulating forward translocation [30]. 

Based on mutagenesis of conserved residues in Sc RNAP 
II, we propose that important aspects of this hinging and 
switching mechanism are strongly conserved in evolution. 
Tt RNAP appears to have two specialized hinges HI and 
H2 within helical segments surrounding the mobile trigger 
loop that are strongly conserved in other RNAPs. Muta- 
tion of the HI and H2 trigger loop hinges in Sc RNAP II 
strongly supports their importance in transcription. Most 
HI substitutions are loss of function mutations, but substi- 
tutions in and around H2 tend to be gain of function 
mutations that can be rapid in elongation and also error- 
prone (Figure 5). Hinges H3 and H4 clustered toward the 
N-terminal end of the bridge helix appear to induce bend- 
ing modes Bl (straighter) and B2 (more strongly bent) to 
support catalysis and translocation. Essentially all RNAPs 
appear to have a flexible H3/H4 bridge helix hinge region, 
but H3/H4 glycines may be differently distributed in differ- 
ent RNAPs, indicating that somewhat different bends, 
flexibility and dynamics may be emphasized for different 
RNAPs. Because bridge helix flexibility and bending 
appears to contribute to multiple steps of the bond 
addition cycle, indicating multiple bending modes, full 
conservation of H3/H4 across RNAPs of different classes 
and species is not expected. The wide distribution of ap- 
parent switch residues in Tt RNAP appear to indicate a 
long range conformational change associated with trigger 
loop closing and opening and progression through each 
bond addition, pyrophosphate release and translocation 
cycle. 

Mutagenesis of a conserved switch that may have a role 
in pyrophosphate release, located about 20 A from the 
trigger loop HI hinge, supports the importance of switch 



residues in transcription elongation. Some other switch 
residues that are predicted from simulations are specific 
to bacterial RNAPs, and those that have been altered ap- 
pear important for transcriptional functions. In some 
cases, the severity of RNAP mutations appears to reflect 
the roles of residues predicted from our simulations. As 
an example, Ec RNAP |3' K781A, corresponding to key 
switch residue Tt |3' K1079 (Figure 4), is strongly defective 
for transcription compared to surrounding substitutions 
(Table 1) [35]. For multi-subunit RNAPs, therefore, mo- 
lecular dynamics simulation can be used to predict con- 
formational changes in the RNAP TEC and the potential 
functions of key amino acid residues. 

In simulation of the closed-Mg Tt RNAP TEC, we 
constructed extra bonds between the Mg atoms, NTP 
phosphate oxygens, the 3'-0 of the RNA chain and 
catalytic aspartic acid residues. These modifications 
appeared to tighten trigger loop closing and to support 
reasonable catalytic geometry (Additional file 1: Figures 
S1-S2, Movies M1-M2). Another group sought in simu- 
lations to tighten Sc RNAP II trigger loop closing by 
protonation of a trigger loop histidine (Sc RNAP II Rpbl 
H1085; Tt RNAP |3' H1242) [25,26], but we did not 
adopt this strategy because we were unsure how to pre- 
dict the appropriate protonation state of this histidine in 
closed and open RNAP TECs. During simulation, bond- 
ing Mg atoms and enhancing trigger loop closing 
resulted in strongly restrained forward translocation and 
excluded water from the RNAP active site (Additional 
file 1: Figures S18-S20). It appears that our closed-Mg 
simulation provides a reasonable initial model for a cata- 
lytic Tt RNAP TEC and an appropriate comparison to 
closed and open RNAP TEC simulations. 

The simulation of an open trigger loop TEC provides 
a comparison to the closed-Mg model but was not set 
up as a perfect representation of a translocating TEC. 
For one thing, an ATP substrate is accurately loaded in 
the active site, which may restrain forward translocation. 
Normally, a translocating RNAP TEC would be a prod- 
uct complex with an open trigger loop, an incorporated 
NMP and a released pyrophosphate (see Figure 1). Com- 
paring the simulations, the open trigger loop TEC has 
more relaxed trigger loop hinges HI and H2 compared 
to the closed-Mg TEC, for which much more action of 
the HI and H2 hinges is evident (compare Movies Ml 
and M2). Interestingly, a sharp bridge helix H3/H4 B2 
bend develops very early during simulation of the open 
trigger loop RNAP TEC and is stably maintained with- 
out observed B2< >B1 or B2< > initial conformation 

(205I or 2PPB) oscillations. The B2 bending mode is 
expected to generate mechanical translocation force on 
the RNA/DNA hybrid, which appears to be indicated by 
strong forward thrust of the DNA/DNA duplex and ex- 
pansion of the RNA/DNA hybrid during open TEC 
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Table 1 Summary of mutations. indicates sensitivity to 6-azauracil 



Structural location 


Subunit, residue and mutation 


Activity in Sc (Tt or Ec) 


References 




T+ DM AD 

It KINAr 


Cr- DM AD 

tC KINAr 


C/- DMAD II 

bC KINAr II 


in vivo 


in vitro 




Switch residue 


p njj/ 


p nO/oL. 


Rnh9 R7^^n/A 
npuz n/ODvj/A 


LtrLl lai 


[CL IIIIIIUl (JclcCLbj 




Switch residue 


P N683 


P N811 


Rpb2 N834A 


6-AU' 


Very slow, moderate fidelity defect 


this paper 


Mg^^ binding 


P E685 


P E813 


Rpb2 E836A 


6-AU' 


Slow, severe fidelity defect 

VCL olUVV clUI ly d LIUI 


this paper 

L4U,4z,4o,4yj 


2+ 

Mg binding 


p uoou 


p uo \ ^ 


Rnh7 nR^7A 

npuz VJOD / r\ 


1 Pthpil 
LcLI Idl 


(^r clr^\A/ olr^nn 3"l"ir^n^ 
VUL blUVV clUI ly d LIUI 1^ 


in IS paper [4z,4oj 


Switch residue 


P N872 


P N1099 


Rpb2 1013A 




Slow 


this paper 


Y-phosphate 


P S878 


P SI 105 


Rpb2 S1019F 


6-AU' 


Slow, severe fidelity defect 


this paper 


Y-phosphate 


P R879 


P R1106 


Rpb2 R1020K/Q 


Lethal 




[43] 


Switch residue 


P' D784 


P' D505 


Rpbl D526A 


Lethal 




this paper 


H3/H4 hinges, bridge helix 






Rpbl G819A, 


Lethal 




this paper 


rij iiiiiyc, uiiuyc iiciix 


R' Gl 076 


P' G778 


Rnhl GRIQA 




"^inw minor firlpiitx/ Hpfprt 


Llllb paptrl , L 1 ^,3JJ 


nj lllliytr, UllUytr IIcIIa 


P' Al 077 


j3' A779 


Rnhl GR70A 




'^I^\A/ rnnHpr;^tp fiHplit\/ Hpfprt 


thic nanpr flA^'^l 
Llllb pdptrl, L 1 i;^-'J 


rij/n^ iiiiiycb, uiiuyc iiciix 


R' Rl 07R 

p n 1 u/ o 


R' R7R0r 
p n/ ou^ 


Rnhl RR7in/A 




RR7in c:|n\A/ RR71A lpth?il 
noz. 1 oiuvv, noz. i r\ icli idi 

(Ec termination defects) 


Llllb pdptrl, IDUj 

[14,35] 


H3/H4 hinges, bridge helix 


P' K1079 


P' K781A 


Rpbl E822 




(Ec very slow elongation) 


[35] 


H4 hinge, bridge helix 




P' G782 


Rpbl G823A 




1 ittlp pffpff rrinHpratp fiHplitv Hpfprt 


this paper, [14,35] 


n 1 iiiiiyc, Liiyytri luup 






Rpbl G1073A, P1075A 


Lethal 




thic PiapiPr 
Llllb pdptrl 


HI hinge, trigger loop 


P'G1230 


P' G924 


Rpbl G1073A 




Slow, strong fidelity defect 


this paper, [2] 


HI hinge, trigger loop 


P' E1231 


P' E925 


Rpbl E1074A 




Fast 


this paper, [2] 


HI hinge, trigger loop 


P' PI 232 


P' P926 


Rpbl P1075A 


Lethal 




this paper, [2] 


HI hinge, trigger loop 


P' T1234 


P' T928 


Rpbl T1077A 


Lethal 




this paper, [2] 


HI hinge, trigger loop 


P'Q1235A 


P' Q929 


Rpbl Q1078A 


Lethal 


(Slow, error-prone in Tt) 


this paper, [2,29,33] 


H2 hinge, trigger loop 






Rpbl G1097A, P1099A 




Rapid in elongation, moderate 
fidelity defect 


this paper, [2] 



simulation (Figures 2 and Additional file 1: S3-S4). Based 
on the enhanced bending of the bridge helix (B2 mode) 
observed for the Tt RNAP open trigger loop TEC com- 
pared to a straighter bridge helix in the closed TEC (Bl 
mode) (Figure 4, Additional file 2: Movies Ml and Add- 
itional file 3: M2), we support the proposal previously 
made by others that bridge helix bending against the 
RNA/DNA hybrid may help to stimulate forward trans- 
location [12]. Because the translocation step is short 
(about 3.4 A), limited motions of the bridge helix followed 
by nucleic acid sliding may largely account for the precise 
forward step of RNAP. The proposed translocation mech- 
anism should limit RNAP to single base steps because the 
mechanical contribution to the step appears to be frac- 
tional, making hyper-translocation unlikely. Nucleic acid 
sliding in the presence or absence of the incoming NTP 
could account for the remainder of the forward step. If the 
incoming NTP is present, the trigger loop is expected to 
close, preventing hyper-translocation. 

Based on simulations, it appears that a set of hinges HI 
to H4 within helical segments surrounding the trigger loop 
(HI and H2) and proximal to the N- terminal end of the 



bridge helix (H3/H4) are important in catalysis and trans- 
location. Somewhat surprisingly, combining mutations in 
multiple RNAP hinges appears to give additive rather than 
synergistic transcriptional effects (Figure 5), although exo- 
nuclease III mapping may indicate somewhat greater func- 
tional coupling of hinges, at least in the presence of an 
incoming NTP (Figure 6). So far, functional assays appear 
to indicate that hinges H1-H4 function somewhat inde- 
pendently during the elongation bond addition cycle. 
Hinges H3 and H4 are very close together on the bridge 
helix, and it is not yet clear whether H3/H4 should be con- 
sidered to be one complex hinge or two semi-independent 
hinges. In the future, other bending modes of the bridge 
helix may be detected as new simulations and x-ray crystal 
structures become available. 

Conclusions 

Using all atom molecular dynamics simulations, compari- 
son of Tt RNAP TEC structures was used to generate a 
model for cycling between catalytic and translocating con- 
formations during elongation. The catalytic TEC appears 
to have a tightly closed trigger loop and, as a result of loop 
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closure, a relatively dehydrated active site. Tight trigger 
loop closure appears to be facilitated by GXP hinges 
within a-helices at the base of the loop. Near the N- 
terminal end of the bridge a-helix, a complex glycine 
hinge (s) also appears to support the catalytic RNAP struc- 
ture by maintaining a relatively straight but somewhat 
extended bridge helix conformation. 

By contrast, the open conformation of the trigger loop 
supports translocation. In the open conformation, GXP 
trigger loop hinges are relatively relaxed, and the more C- 
terminal GXP hinge (H2) becomes the site of helix break- 
ing, increasing the size of the loop. During simulation of 
the open trigger loop RNAP TEC, the bridge helix N- 
terminal hinge (s) bends against the RNA/DNA hybrid 
stimulating forward translocation. In simulation, switching 
between catalytic and translocating RNAP TECs is sup- 
ported by amino acid residues that trade atomic contacts. 
Mutation of hinge and switch amino acid residues and 
assay of mutant RNAPs supports the proposed model. A 
conserved switch is identified by simulation, mutation and 
assay that appears to function in pyrophosphate release 
after catalysis and trigger loop opening. 

Methods 

Tt RNAP TEC simulations 

Molecular dynamics simulation of Tt RNAP closed trigger 
loop TECs (PDB 205 J with some missing chain taken 
from PDB 205I) in explicit solvent with counterions was 
done as previously described [13]. AMPcPP (a-, |3- 
methylene ATP) in the original crystal structures was sub- 
stituted with ATP for simulation. During modeling, Tt 
RNAP |3' H1242 was not protonated because we were not 
certain how to predict the protonation state of this residue 
[25,26]. Thio-substitutions in nucleic acids used to block 
RNA cleavage by Tt RNAP in crystals were not present in 
the protein data bank files and were not considered in 
modeling [16,21]. For the closed-Mg simulation, extra 
bonds were assigned according to reference [16]. For this 
simulation, Mg-I was covalently bonded to the RNA 3'0, 
|3' D739 ODl, |3' D741 ODl, p' D743 OD2 and the a- 
phosphate O. Mg-II was covalently bonded to the a, |3 and 
y-phosphate oxygens and to D739 02. Bond lengths 
were taken from the original 205J crystal structure and 
were maintained with a force constant of 200 kcal/mol/ 
A^. For the closed TEC simulation, extra bonds to Mg 
atoms were omitted. Simulation of the open trigger loop 
TEC (PDB 2PPB) was done by essentially the same proto- 
col as the closed trigger loop TEC. For consistency, strep- 
tolydigin, which was present in 2PPB but absent from 
205J, was omitted in simulation of the open TEC. Missing 
chain in the 205J closed trigger loop TEC cleft region was 
patched with chain from the 205I TEC structure using 
the program MODELLER [51]. Simulations were done on 
the NSF XSEDE (Extreme Science and Engineering 



Discovery Environment) national computer system using 
the University of Texas, Austin, Sun Constellation Linux 
Cluster TACC (Texas Advanced Computer Center) Ranger 
computer. To parallelize large computer jobs, simulations 
were done using the program NAMD (Not (just) Another 
Molecular Dynamics Program) [24] using AMBER 
(Assisted Model Building with Energy Refinement) [23] 
force fields. To limit computation time, a and co subunits 
were omitted from the simulated files, as previously 
described [13]. Because no large distortions were observed 
near omitted a and co subunits, these deletions in the TECs 
did not appear to negatively affect simulations. Potentially, 
however, conformational changes observed in simulations 
may be accelerated or amplified by omission of a and co 
subunits. Additional simulations with all RNAP subunits 
have now been run but for simplicity of presentation are 
not reported here. These simulations produce closed-Mg, 
closed and open RNAP TEC conformations that are very 
consistent with those shown here. 

Equilibration consists of 500-1000 steps of minimization 
followed by instantaneous heating to 320 K and a gradual 
release of harmonic restraints from 5 to zero kcal/mol/ A^ 
over 500 ns. This is followed by 10.5 ns of production 
simulation run at 2 fs/step with the SETTLE algorithm and 
TIP3P water. A Langevin thermostat and barostat are used 
to maintain 320 K temperature and 1 atm pressure. Simu- 
lations were run using NAMD version 2.8 and with Amber 
fflO force fields from Amber 11. Non-bonded (electrostatic 
and Van der Waals) forces are calculated out to 10 A with 
switching at 8 A for smoothing the cutoff. For the calcula- 
tion of electrostatic forces. Particle Mesh Ewald summation 
is used with periodic boundary conditions [23,24,52]. 

Molecular images were prepared using the program 
Visual Molecular Dynamics (VMD) [53]. Representative 
PDB files from simulations will be made available on re- 
quest to ZFB (burton@cns.msu.edu). 

Strains and in vitro mutagenesis 

Site-directed mutagenesis, production and purification 
of Sc RNAP II mutant proteins was done essentially as 
previously described for rpb2 mutants [43]. For produc- 
tion of rpbl mutant strains a similar set of plasmids and 
strains was constructed. Phenotypic testing of yeast 
strains was as described [43]. 

In vitro transcription 

Bulk elongation assays were done as previously described 
[54]. Briefly, TECs were assembled by annealing - RNA 
oligonucleotide (5' AUCGAGAGG 3') to the template 
DNA oligonucleotide (5' GGT TTG CCC CGT TGG 
ACG TGT GGA TTG GGA TGG CTA TTC GCC GTG 
TCC CTC TCG ATG GCT GTA ACT ATC CTA TAG 
G3'), and incubation of Sc RNAP II with the RNA-DNA 
hybrid followed by addition of the non-template DNA 
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strand (5' CCT ATA GGA TAG TTA GAG GGA TCG 
AGA GGG ACA GGG GGA ATA GGG ATC GGA ATG 
GAG AGG TGG AAG GGG GGA AAG GGT A 3 ). The 
RNA in the resulting TEG was labeled by 2 min incubation 
with 0.15 \iM a-^^P GTP, and the resulting GIO TEG was 
purified from the unincorporated GTP and nucleic acids 
using the centrifugal filtration through a 100 kDa cutoff 
membrane. Transcription was done with 200 [iM NTPs, 
and the time points were 5, 10, 20 and 40 s. For most 
mutants the experiment was done at least twice. Error bars 
represent standard error from a single phase exponential 
fit of the runoff accumulation time course. 

Fidelity assays 

Gompetition fidelity assays were done as in Walmacq 
et al, 2012 [32] and as shown in Figures S9 and SIO. 
Accurate GMP incorporation at GIO is discriminated 
from misincorporation of AMP at A* 10 to indicate the 
misincorporation product) by the different mobilities of 
the short RNAs on gels. The results of this assay and the 
calculated fidelities are very similar to those obtained 
using more complex normalization protocols [20] (Add- 
itional file 1: Figure SIO). 

Exonuclease III footprinting of RNAP translocation state 

Sc RNAP II was isolated and assembled into G9 TEGs, as 
described [20,43]. TEGs were extended to GIO and All 
with addition of 10 [iM GTP and 100 [iM chain terminator 
3'-dATP. Ec exonuclease III (New England Biolabs; 3.3 U/ 
|il) mapping of the upstream Sc RNAP II TEG border was 
done as described in Kireeva et al. [20]. Exonuclease III 
digestions were done at 25 G. DNA and RNA were resolved 
in 10% polyacrylamide-7 M urea gels. Quantification of ^^P 
in wet gels was done using a Typhoon 9200 (Amersham 
Biosciences) phosphorimager and ImageQuant software. 
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Additional files 



Additional file 1: Supplementary figures and legends. 

Additional file 2: Movie Ml. Closed-Mg (closed trigger loop; bonded 
Mg atoms) Tt RNAP TEC {Thermus thermophilus RNA polymerase ternary 
elongation complex). DNA is gold. RNA is red. ATP is red (stick). The 
bridge helix is green transparent. The trigger loop is green opaque. 
G1230 (HI), G1233 (HI), G1255 (H2), G1076 (H3), G1080 (H4) and G1081 
(H4) are shown in space-filling representation. Bridge helix bending 
mode Bl is established and held stably through the simulation. At the 
beginning of the movies, closed-Mg and open RNAP TECs are very 
similar in conformation except for the trigger loop conformation. 

Additional file 3: Movie M2. Qpen (open trigger loop) Tt RNAP TEC. 
DNA is gold. RNA is red. ATP is red. The bridge helix is blue transparent. 
The trigger loop is blue opaque. G1230 (HI), G1233 (HI), G1255 (H2), 
G1076 (H3), G1080 (H4) and G1081 (H4) are shown in space-filling 
representation. Bridge helix bending mode B2 is established rapidly and 
held stably through the simulation. At the beginning of the movies, 
closed-Mg and open RNAP TECs are very similar in conformation except 
for the trigger loop conformation. 
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